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to compute the effective permittivity tensor (in the quasi-static limit) of this composite in a static 
magnetic field. In order to find the scaling exponents near the percolation threshold Pc, we use a 
high contrast or low frequency expansion combined with scaling. The results of the two methods 
are in agreement and predict that near pc (and below, that is, in the dielectric region), the Faraday 
effect is greatly enhanced. 
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I. INTRODUCTION 

We will concentrate on the problem of the effective permittivity tensor of a binary mixture in a static magnetic field. 
Since the Faraday effect is usually very weak in dielectrics, we study a mixture made of a dielectric host matrix (with 
a negligible Faraday effect) with metallic inclusions. In such a system, the Hall effect in the inclusions is expected to 
induce a measurable Faraday effect in the composite dielectric even when this effect is totally negligible in the pure 
dielectric host. We will study the quasi-static limit, the wavelength of the incoming wave being much larger than 
the typical inhomogeneity length (which is for example the size of the metallic grains). Therefore, the composite 
medium can be seen as quasi-homogeneous, and the equivalent homogeneous material is called the effective medium. 
Determining the effective medium properties of disordered materials (such as composites, suspensions) is a difficult 
problem and one has often to resort to perturbative methods (low field, low density or low contrast) which cannot 
be applied for high magnetic field for instance. We propose here two different approaches to this problem. The 
first one is based on the use of the replica method together with a variational principle. This treatment possesses 
the advantage that it is non-perturbative and may thus be useful for strong disorder ror strong fields. This method 
was successfully applied to different problems (the random resistor network problemlil. Hall effect in compositesH). 
Moreover it has been shown that it can give reasonable values for the permittivity tensor if the system is not too close 
to the percolation threshold pc- However, the critical exponents are not reproduced correctly (one usually gets mean- 
^ ' field or effective- medium-approximation exponents). In order to present an alternative discussion of this problem, 
, and obtain the correct exponents near pc, we introduce a high contrast expansion. This is essentially an expansion in 
' powers of the ratio of resistivities or permittivities of the two components, which can be made very small by making 
5^ , the frequency of the incoming wave very small. This expansion can be used for weak magnetic fields as well as for 
strong magnetic fields. In order to discuss the critical properties near pc, we iSjPjply some scaling ansatzes to that 
expansion, which are based upon previous discussions of d.c. magneto-transportau. 

The Hall effect in percolating composites has beierL.s.tj|idied using a number of different methods which enabled 
the critical behavior to be determined quite reliablyETtrffl. However, until now the Faraday effect in metal-dielectric 
composites was only discussed using a Clausius-Mossotti-type approximatioaEl, which is good for dilute systems, and 
a Bruggeman-type self consistent effective medium approximation (SEMA)lj, which exhibits a percolation threshold 
but with incorrect values of the critical exponents. In the present study we employ different approaches (see above). 
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Both of our approaches are not hmited to dilute systems, and one (the high contrast expansion together with scaling 
ansatzes) is expected to lead to reliable results for the critical behavior near pc- 

Let us first recall some facts about the Faraday effect. When an isotropic material is subjected to a static and 
uniform magnetic field B directed along the z-axis, it can be described by the following permittivity tensor: 

(1.1) 

where e must depend upon B (see, e.g., Ref. |lo| ). In a homogeneous medium, the dispersion equation gives rise to 
two solutions which are the left and right circularly polarized waves with two different refractive indices (e and e are 
positive and real and e > e, which ensures that the wave is undamped) 




n± = \/e±e. (1.2) 

If a linearly polarized wave of frequency lu propagates over a distance L through this medium, the polarization plane 
will rotate (the so-called Faraday effect) by an angle (c is the light speed in vacuum) 

e = -L{n+-n_). (1.3) 
c 

Usually e is very small compared to e (i.e., the Faraday effect is weak) and the rotation angle is therefore approximately 
given by 

e^^-L± (1.4) 

c V£ 

We will also use the Faraday coefficient, which is defined hy T = e/^/e. We note here that the Faraday effect is 
usually weak {T ranges from 10~^ in dielectrics to 10^^ for thin films of metallic iron). 

We will study the case where the material is a random binary composite medium made of a dielectric host with 
metallic inclusions, and where a static uniform magnetic field B is applied along the z-axis. We suppose that the 
medium has a position dependent permittivity tensor e(r) which is an independent random variable at each point r 
(uncorrelated at different positions), distributed according to the following probability density 

p{e) = pS{i - sm) + (1 - p)S{i - ii). (1.5) 

Let us note that in real materials, the grains have finite sizes and that in a finite frequency calculation one should take 
this into account. Ho wev er, in the quasi-static limit, the grain sizes are irrelevant and this simplified characterization 
of the disorder [Eq. (pT^)] is justified. 

Equivalently, the local resistivity tensor p{r) (related to the permittivity tensor e hy p — 4^) is a step function 
that is equal to pM inside the metal and to pi inside the dielectric component. It will be represented, with the help 
of the appropriate characteristic functions Oui'r), ^/(r), as 

/5(r) =pA/^M(r)+/5/0/(r), (1.6) 

„ , ^ /, / X f 1 for r inside the metal, ,^ _s 

0M(r) = l-^?Kr) = |o otherwise. (^-^^ 

In the dielectric component, the permittivity is taken to be 

ii ^ eii. (1.8) 

where e/, the dielectric constant of the host, is a real scalar quantity and is independent of B. The metallic component 
is non-percolating, and is characterized by a free-electron-like resistivity tensor in the presence of a magnetic field 
B II z or, equivalently, by its permittivity tensor 




eM=£M\ -TTTp- -rpw , = cJcT cx |B|, (1.9) 



where uJc is the cyclotron frequency, r is the conductivity relaxation time, and Em = is purely imaginary and 

independent of B ((Tm is the conductivity of the metallic component). We assume that this form continues to be valid 



2 



even at finite frequencies. This probably means that the entire subsequent discussion will not be valid for optical 
frequencies in the visible range. But it will be relevant for frequencies up to, and including, the microwave regime. 

These assumptions mean that the host exhibits no intrinsic Faraday effect and the metallic component has no 
intrinsic magneto- resistance, only a Hall effect. In terms of resistivity, the metallic component is characterized by a 
free-electron-like resistivity tensor (obtained by inverting em) 




PM = Pm\ -H 1 , (1.10) 

and the impedance of the dielectric host is an imaginary scalar tensor 

47r 

pi = pil, pi = . (1-11) 

lUEi 

We will assume that the quasi-static approximation can be used (i.e., both wavelength and skin depth are large 
compared to the sizes of metallic inclusions) . 

We define the bulk effective permittivity tensor ie of the medium by the following relation 

(e(r)E(r))=4(E(r)), (1.12) 

wh ere E is the electric field and where the brackets denote a quenched average over the probability distribution given 
by ( |l.5| ) (or equivalently a spatial average over the volume of the sample) . The effective medium will be homogeneous 
and isotropic and we expect an effective permittivity tensor of the form 

(1.13) 

We can also define the bulk effective complex resistivity tensor by 

Pe-(J) = (pJ), (1.14) 

where J is the local current density. 

In order to evaluate the effective properties of the heterogeneous medium, we can proceed in different ways. The 
first one (which will be presented in Section ||) relics on the observation that the effective permittivity tensor can be 
(exactly) related to the inverse of a random operator M. The problem is thus reduced to the calculation of (M~^), 
and we will use replicas (and a variation al ap proximation) to ev aluate this quantity. Let us note here that it is in 
principle equivalent to compute Eg using ( 1.12| ) or using (1.14), since the product of these tensors is proportional 




to the unit tensor /. However, since we use approximations, the two procedures are not necessarily equivalent. In 
fact, it has been shownEj'EHI that for the variational approximation, computing or leads to different results, and 
that it is more reliable to calculate the permittivity (or the conductivity) than it is to calculate the resistivity. 

Since the variational method is non-perturbative, it possesses the advantage that it is reliable in the whole range of 
concentration and for any strength of the disorder. However, near the percolation threshold p^, this method leads to 
mean- field exponents which are usually not accurate (especially in three dimensions d — i). Moreover, it is difficult 
to control the quality of the variational approximation. 

In order to describe the Faraday effect near pc, we will use a second approach, which relies on the fact that if w is 
small enough, we have two small parameters puj Pi and Hpm/pi- The local electric field E(r) and current density 
J(r) can be found by defining a vector potential A(r) such that 

J(r) = VxA(r), (1.15) 

and then solving the equation 

VxE = (1.16) 

using the constitutive relation E(r) — p(r)J(r). It is of course out of the question to solve this equation exactly, 
therefore we will expand A in powers of the contrast between the two components, i.e., the resistivity ratio puj pi- 
This approach is valid only in the quasi-static regime. Besides confirming results obtained by means of the replicas, 
when combined with scaling it allows us to predict the behavior of the medium near the percolation threshold. 

The rest of this paper is organized as follows. In Section || we apply the replica method to the Faraday effect in 
random composites. Results are obtained for the scaling behaviors near the percolation threshold. In Section 111 we 
present the high contrast or low frequency expansion. Section presents a scaling theory based upon the approaches 
described in the two previous sections. Section ^ summarizes the main conclusions from this work. 
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II. REPLICA APPROACH 



In this section we present the rephea approach. We compute here the effective permittivity tensor of a binary 
mixture, where the tensor e is a random variable equal to ei with probability p and to £2 with probability q — 1 ~ p 



(each component haying a permittivity of the form given in equation (1.1)). The calculation is essentially the same 
as for the Hall effectfi and we give the main steps of the derivation in Appendix A. We obtain the following equations: 
The first one gives the longitudinal effective permittivity 



d^e-" ^ V . (2.1) 







where the brackets still denote an average over ( |1.5| ), and we also get two coupled equations for the transverse 
permittivities £e, £e'- 

q roo 

I = -7: due-"{ln(e-"^/3^=) +ln(e-"^/3^^)}, (2.2) 
2 Jo 



due~"{ln( 



ln( 



)}, 



(2.3) 



where X = e — e and Y = e + e. Eq. (2.2) determines Eg self-consistently while (2.2) determines ie after Se is known. 

Eq. {2A_) is the same as that obtained in Refs. |^ and |l^ for the bulk effective conductivity of a binary mixture 
of zero-field conductivities azi with concentration p and az2 with concentration 1 — p. That equation was studied 
thoroughly in Ref. [ij; it displays a percolation t hre shold at Pc = 1 ^ e^^/^ ~ 0.28. 

The Faraday effect is contained in ( pT^ ) and (2.3). We first note that the percolation threshold is independant of 
the magnetic field. Let us note here that the percolation threshold, which is a geometrical quantity, is still meaningful 
here since we are working in the quasi-static limit. 



By expanding the logarithms (these expansions are valid for p < 1/2), we obtain from (2.2) 



1 



3ee 



1 



3ee 



(2.4) 



and from (2.3) 



Ee = £2 + {£1 - £2) A"(- 



\n+l 



n>l 



[1 + ^{X,-X2)][l + ^{Y,-Y,)r 



(2.5) 



where A = p/q, with p the fraction of component 1 and q = I — p the fraction of component 2. By some simple 
algebraic manipulations, these equations can be rewritten as 



o = \+inq-^ + Y.-i-r^' 



n>l 



(1 



1 



ra(ai— 02) 
3ae 



71(01—02) ^2 



30e 



)^ + C 



3ae 



(2.6) 



and 



/?e = /?2 + Wl 



n>l 



n+1 



(2.7) 



where and f3i describe, respectively, the diagonal and the off-diagonal transverse elements of the conductivity tensor 
of component i [ai = (Ji/{1 + Hf), Pi — aiHi/{l -\- )], and where <Ji is the zero-field conductivity of component i 
[Hi is the Hall-to-Ohmic resistivity ratio in component i: it is proportional to the magnetic field strength |B| — see 
(L9)]. The quantities ae and (ie are the effective coefficients of the composite. Eqs. (^) and (2/7) are identical to 
the equations obtained in the case of magneto-transporto. We have thus shown that the results obtained in that case 
can be continued analytically to complex values of a and (3 (we can go from the Hall effect to the Faraday effect by 
changing a to £ and /? into ie). 
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We now study these equations in the case of a metal-dielectric mixture (the metal will be component 1 and the 
dielectric component 2). We assume here that the Hall effect in the dielectric is very weak so that a-ijoLi <C 1 or 
Hi ~ 0, and we use H instead of H\. We also assume that cte ^ ai. 

We will first study the weak field regime H <^ 1, in this case Pi ai. Near the percolation threshold Pc — I — e~^l^ 
of component 1 (Ap = p — Pc ^ 1), we find that ae has the following scaling behavior a,, ~ ai|Ap|0( °^p°^ ), where 
the scaling function 0(z) satisfies 

P^_Ap_ a2 +\Ap\M, (2.8) 
where A = 3 J2n>i ^(~)"^^ (with Ac = Pc/l — Pc)- We thus obtain the following equation for (j) 

A0^-^-|=O, (2.9) 
where e = +1 if p > and e = — 1 for p < pc, and where z = ^^^p^ ■ The solution of this equation is <f) = 



jxi-^ + J ^ + ^4^) I which for small z becomes (up to a constant factor qd'i) 



(z) cx <j 3^^^^' P^P" (2.10) 
^ ' I Const; p> Pc- 



We can now easily obtain the behavior of /3e from (EJ) (for p < pc) 



/3e OC ^Qfg CX , (2.11) 



which is proportional to H. 

We now con side r the regime H ^ 1 for p below pc (which is the interesting regime for the following) . In this regime 



Pi ^ a,;, and ( |2.6| ) then reads 



0- --T^ + ^A 2.12 

qc 3ae Pf 

from which we can deduce that the scaling behavior of Ue is of the form ae — ^\Ap\(l){-^^^), where the scaling 



fmiction 0(z) has a behavior like that of (f){z), up to the same constant factor qc/3 [see (2.10) 



The behavior of P^ is a priori different, since now Pi oc jj- ^ ai (x The equation for P^ then becomes 

A" 1 

^ /5i ^ -i(-)"+ \. ,„ (2.14) 
,T>1 " (/3i/3ae)^ 



which for p < pc leads to the following behavior 



which is again proportional to H. 



We will discuss the physical consequences relevant to the Faraday effect in Section IV 



III. LOW FREQUENCY OR HIGH CONTRAST EXPANSION 

In this section we derive an expansion for in powers of the complex resistivity ratio puj Pi- Throughout this 
section we will assume that to is small enough so that both pi 3> pu and pi ^ Hpu-, and that the quasi-static limit 
is valid. 
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Under these assumptions, the local electric field E(r) and current density J(r) can be found by defining a vector 
potential A(r) such that 

J(r)=VxA(r), (3.1) 

and solving the equation 

= VxE = Vx {p{r) ■ [V X A(r)]} (3.2) 

along with appropriate boundary conditions on n x A(r) at the system surface (n is the unit normal vector to that 
surface) . 

We recall that the local resistivity tensor p{r) is a step function equal to pM inside the metal and to pi inside the 
diele ctric component, and that it can be represented with the help of the appropriate characteristic functions as in 

In connection with ( [3.2[ ) it is useful to define a Green tensor G^f^r, r') by the following equations 



{V X [p • (V X G^;;0]}a - k^Gl'> = 6o^pS^{r - r') for any r, r', (3.3) 

n X G^p"^ = for r at the system surface. (3-4) 

This tensor can be used to solve ( ^.2[ ): 

A(r) = A(")(r) - / dr' lim[V' x G('')(r,r')] • /5(r') • [V x A(*')(r')], (3.5) 

where A'^*'^(r) is a vector field that satisfies the same boundary conditions as A(r), but is otherwise arbitrary. Note 

that we need to use the limit fc of V' x &P\r, r') here: We could not take that limit in (3.3), because then the 
equations for would have no solution (see Appendix B for a discussion of this point). 
Since we intend to expand A in powers of pAi/pi, we define 

G^'^ EE lim G^P\ (3.6) 



It is then possible to transform (3.3), (3.4) into an integro-differential equation that relates and 

Gi%v') = gI^IM-J dr"[V" X Gl^Hry')]-PM0Mir") ■ [V" x G^;\v",v')]. (3.7) 

Iteration of this equation leads, in the usual way, to an expansion of in powers of pM around We note 

that, although G'"^^{r, r') is not a symmetric kernel (because p(r) is a non-symmetric tensor), G|/^(r, r') is symmetric 
because pi is symmetric (in fact, pi is a scalar tensor) [see Appendix B for a discussion of this point] 

GS(r,r') = G(^j(r',r). (3.8) 

A possible choice of A'^'^^ in (^]^) is limpJ^,,^o A — henceforth we adopt that choice. If we then take the limit pM 
also in the functions A(r), &f\r,r'), we conclude that 

dr' lim [V' X G(^) (r, r')] • Pi9iir') ■ [V x A^") (r')] = 0. (3.9) 

k — '0 

Using the above mentioned power series expansion for &p\ this result can be extended to hold also when G^^' is 
replaced by : 

I dr' lim [V X G^") (r, r')] • p/^/(r') • [V x A^") (r')] = 0. (3.10) 

J k^o 

Using this result together with (3.7), (^.5|) can be transformed into an expansion for A(r) in powers of pM- 
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A(r) = A('')(r) - / dr' lim[V' x G('''(r,r')] • PM^Mir') ■ [V x A(")(r')] 

= A("' (r) - f dr' lim [V x G^^'> (r, r')] • PM^Mir') ■ [V x A^o) (r')] 
+ 0(pL). 

We also note that A(°) (r) can be obtained by an expression that is the analogue of ( |3.5[) , namely 
A(o)(r) = A("°)(r) - / dr' lim[V' x G(^)(r,r')] • /5/6'/(r') • [V x A("°)(r')]. 

J *:^0 



(3.11) 



(3.12) 



This was obtained by replacing A'''^(r) by A(°'''(r) in (3.5), and then taking the limit Pm ^ in that equation. Eq. 
(3.12) is especially useful if we assume the following boundary condition for A(r) and A'^"^(r) 



A(r) — A^''-'(r) = -(e x r) at the system surface, 



and choose 



A^''°''(r) = -(e X r) everywhere, 
3 

where e is some unit vector. This last choice corresponds to a uniform current density 



(3.13) 



(3.14) 



V X 



-(ex r) 
3^ ' 



= e, 



which is equal to the volume averaged current density for both A(®)(r) and A('^'')(r), which satisfy ( 3.13 ), 

(V X a'*')) = (V X A^O'')) = e. 

Wc thus find that A^^''-' is given by 

A("'')(r) = i{exr)- J dr' limJV x G(^)(r, r')] • p/6ij(r') • e. 



(3.15) 



(3.16) 



We recall that the bulk effective complex resistivity tensor of the system is defined by [see (1.14) and 3.1)] 

Pe-(VxA) = (p.(VxA)), (3.17) 



for any A of the form (3^). Using (1^) and ( 3.11 ), we can expand an arbitrary element of the tensor pe in powers 
of Pm (f, e are arbitrary unit vectors) 

(f.pe-e)=f.(/5.(VxA(-))) 

= f • /}/ • (0,(V X A("-))) + f • PM ■ {OMi^ X A(°-))) 

-f -iy" drp,0,(r)-Vx 

dr' lim[V' X G(^'(r,r')] • pM0M{r') ■ [V x A(°")(i-')] 
0{plt). (3.18) 



(note that (6'/(V x A'^^'^^)) is just the spatial average of V X A^"*^) over the subvolume of the dielectric component, 
while (61a/ (V X AfO'^'jWs the avera ge of the same quantity over the metallic subvolume). The integ ration over r can 
be performed using (3^) and ( 3.16 ), leading to the following result for the double integral of ( |3.1^ ) 



J dr'[V' X A(°f)(r') - f] • pMOM{r') ■ [V x A("'')(r')]. 



(3.19) 



Part of this cancels the second term of ( 3.1^ ) , and we finally get 
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f • Pe(/5/,PM) • e = f • • (0/(V X A^O'^))) 

+ ((?m(V X A("f)) • pM ■ (V X A(o-))> + 0{pli). (3.20) 

If the microstructure is isotropic, then since pi is a scalar tensor, Pe{pi, 0) [the first term on the r.h.s. of ( p.20[ )] is 
also a scalar tensor, and it is clearly independent of pM and hence of H. The vector potentials A'''''^ A^'^^', which 



appear in the second term on the r.h.s. of ( 3.20 ), satisfy different boundary conditions at the system surface [see 
(3.15)]. Inside the metallic subvolume, those potentials can also be viewed as resulting from boundary conditions on 
n X A(°^ at the interface between the two components . Those latter boundary values are entirely determined by the 
microstructure when we impose the requirement that the electric potential must be constant over every connected 
subvolume of the metallic component, but the precise local values of A^''-'(r) inside those subvolumes also depend 
upon the Hall-to-Ohmic resistivity ratio of the metal H . Nevertheless, we now argue that even the second volume 
average which appears in ( ^.2C| ) is independent of H in the two limits i7 <C 1 and H ^ 1. The only H dependence 
in those limits arises from the explicit pM factor in that term. 

In order to prove this, we note that if the resistivity ratio pAi/pi is small enough, then the current distribution 
inside the metallic subvolumes, though different for <C 1 and for if ^ 1, will be saturated in both limits: In the 
weak field limit this is obvious, while in the strong field limit this holds in a percolating system whenever the magnetic, 
field dependent correlation length S^h, which diverges as ff — > oo, is greater than the percolation correlation lengthcl 

Recalling that B || z, and assuming that the microstructure is either isotropic or cubic, we now get that the diagonal 
elements of pe are given by 

P^:liPi,PM)^Pi{0i{y X A(""))„) +pm(^m(V X A(0"))2) + 0(p1,), (3.21) 

while the nonzero off-diagonal elements are 

pi'^iPi, Pm) = -Pi^{pi, Pm) = HpM{eM[{V X A(o-)) X (V X A(o«))],) + 0{pl,). (3.22) 

Recalling also that 

47r 

Pe = — , 

we finally get the following results for ie {<Jm = ^/ Pm)'- 

ei f iLoe] (0Af(V X A(o^))2) 



£e{£l,PM) 






1 




(3.23) 



(6l/(V X A(o^))^) AnaM (6'/(V x A(o^)))2 

AiraM {OiiV x A(o^)))2 

iujejH (6Im[(V x A("=")) x (V x A^"?'))], 
~ AnaM (6l/(V x A("^)))2 

This expression is the main result of this part and we will discuss it in the next section. 

IV. SCALING THEORY 

In this part, we will discuss physical consequences of the results obtained by the different approaches, and we will 
use a scaling theory in order to discuss the behavior of the permittivity below the percolatioa-threshold by comparing 
with scaling theories develop ed e arlier for d.c. magneto-transport in a percolating systenJjQlij. 

As we noted earlier, Eq. ( ^.l[ ) is the same as the result obtained in Refs. |l] and Jll|: It is an equation for the 
bulk effective conductivity of a binary mixture which was studied thoroughly in Rcf. |12|; in particular, it displays a 



percolation threshold at Pc — I — e ^/^ ~ 0.28. The Faraday effect is contained in Eqs. (|2.2|) and (2.3) or, equivalently. 



in (2.6) and (2.7). We first note that the percolation threshold for that effect is independent of the magnetic field 



and, as expected, is the same as the threshold found for the conductivity. 
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We can easily check those results in two limiting cases. First, for zerOpHiagnetic field it is easy to see that one 
recovers Hori's equation for the effective permittivity of a binary mixturellj. Then, for low concentration {p <C 1) 



we find that Se — £/ and ie 



,47ro-flf 



H/{1 + H ), and the rotation angle therefore satisfies 6 (x 1/H for high fields 



and 9 oc H for low fields. This agrees to order 0{p) with the low density expansion results for spherical inclusions as 
obtained, for example, from the Clausius-Mossotti-typc approximation of Ref. |[ Wc note that here 9 is independent 
of the frequency and is always very small. 

We now apply our discussion from Section || to the case of a non-dilute metal-dielectric mixture: For the metal, 
cti = (^i'^'h'^) ^^'^ Pi — ctiH, and for the dielectric (which is assumed to have a negligible Faraday effect) a2 = iujei/An 
and /32 — 0. We will concentrate on the critical region near the percolation threshold pc, where Ap = p — pc is small 
{p is the metal volume fraction). We will now consider Ee and ie separately for the weak field regime {H ^1) and 
the strong field regime {H ^1). 

At low fields, we use the scaling result for a^, and thus find that the effective permittivity {se = ) is given by 



{\Ap\ ''i2^(tm|Ap|3 



(4.1) 



for the regime where ^^"^^2 ^ 1 and p < pc- The real part of thus diverges like l/|Ap| (exponent s equal to its 

2 

SEMA value 1), and the imaginary part is proportional to 



From (2.11), we obtain the following result for 



Be OC 



47rCTM|Ap|2 



H. 



In this regime, the Faraday coefficient will be 



T (X 

In the strong field regime, (3i ~ um/H and ai 
that the effective permittivity is given by 



3/2 „ 

47rcrAf |Ap|3/2 ' 



(4.2) 



(4.3) 



gm/H"^- Using the scaling result for Ue in this regime, we find 



^ — 



Ml 



uje, 



3 \\Ap\ 12tt c7Ai\Ap\^ 



(4.4) 



This scaling behavior is the same as what was found above in the weak field regime. The first term should indeed 
be independent of H, since it corresponds to the universal behavior of the d.c. permittivity near the percolation 
threshold. However, the fact that the imaginary parts of (4.4) and (4.1) are the same to the order shown here is 
accidental, and is probably due to the nature of the approximations used (a similar accident also occurs in the SEMA 
results). Indeed, we can see from (3.23) that these parts depend on the current distribution in the metallic inclusions, 
hence they should be different in the weak and strong field regimes. The important physical conclusion is that, in 

both regimes, Se is independent of H . 

In the strong field regime, we obtain from ( ^.15 ) the surprising result 



£e OC 



47rcrM(Ap)2 



(4.5) 



which is again the same scaling behavior as in the weak field regime. In both regimes, the Faraday coefficient thus 
reads 



J- cx 



3/2 tj 

uje/ H 
47r(TM|Ap|3/2' 



(4.6) 



This means that the Faraday rotation (the angle is proportional to is proportional to the applied magnetic field, 
even for strong fields, i.e., when _ff 3> 1. It is thus clear that we can obtain a large value of the rotation angle using 
such a composite. 

We now consider the results obtained by means of the high contrast or low frequency expansion. The scaling 
behavior of the averages which appear in ( 3.23| ) pfi^n be deduced by comparing (3.2C) to scaling theories previously 
developed for d.c. magneto-transport cocfficicntsfla, and from the property that 
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1 r S percolating cluster; p > Pc 
j(o) = V X A^"^ cx <( r e elsewhere; P > Pc (4.7) 



These considerations lead to 



1 r G anywhere; p < pc 



, , f Ap"*; any p < 1 

(^m(V X A(-))^) oc I ^1^^ . ^^^^ ^ ^ ^ (4.9) 

f An^f; any p < 1 

(^m[(V X A(-)) X (V X A("^))].) oc I ; any p i7 » 1 ■ (^'^O) 

Here ~ H'^" is the magnetic field dependent correlation length, cx Ap^"^ is the percolation correlation length, 
and F{z), G{z) are scaling functions, which tend to nonzero constants when z ^ 1, and to asymptotic forms that 
return the Ap exponents to their H values when z 3> 1 (see Ref. ^) : 

, ^ f const, z ^ 1 , ^ 

F z (X 4.11 

y z - , z > 1 

, ^ ( const, z <C 1 , , 

Gzcx ^ 4.12 
I z - , z > 1. 

The values of the critical exponents which appearpkuthe above scaling expressions, as determined by simulations of 
three dimensional percolating network models, areOLlu 

1/ ^ 0.88, t ^ 2.0, s ^ 0.7, g ^ 0.38, , 

1.^ = 0.5, iff ^6.0, gH = 5.0. ^^'^^^ 

Finally, we get the following results for the Faraday coefficient of a percolating mixture below pc 

— Ap-9-^^/^, i?<l 



coef^H (0m[(V X A(o^)) x (V x A(o«))],) 

OC ^ 

AnaM {9iiV x A(o^))^)3/2 ^ 



3/2 rrl + (gH-g}i,H/u (4.14) 



These results are consistent with (4.6), which was obtained using the replica method, where we expect to find the 
SEMA values s — 1, g = g// = 0. We can also derive a number of physic al c onsequences which follow from both 
approaches [from ( 3.23| ) in the high contrast expansion, and from (4.1) and (4.4) in the replica approach]: 



1. The diagonal part of has an imaginary part that is proportional to Loe'j/aM and is independent of H (up to 
0{pm) h^ the high contrast expansion). This means that there will be some dissipation. 

2. In the same order, has an antisymmetric part which is imaginary and proportional to uie'jH/aM ■ These 
results are valid both for H < 1 (weak field) and for _ff > 1 (strong field), as long as both ujei ^ ctm and 

We note that the scaling predictions for magneto-transport in a percolating system(i|ave been tested experimentally 
only for weak fields, and only in systems that were above the percolation threshold Pc&B. Measurements of the induced 
Faraday effect in a metal-dielectric mixture below pc could therefore provide an important test of those predictions. 



Another prediction which follows from (3.23) is that the induced Faraday effect in a non-conducting metal-dielectric 
composite, which is not necessarily near any percolation threshold, is linear in H = lOcT even when iJ 3> 1, in 
agreement with the replica approach near pc- 
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V. CONCLUSION 



In this paper, we studied the Faraday effect (in the quasi-static hmit) in a composite consisting of a dielectric 
matrix (with neghgible Faraday effect) with metaUic inclusions (which have only Hall effect and no intrinsic magneto- 
resistance). We presented two different approaches leading essentially to the same conclusions. The first approach 
relies on the replica method and allowed us to derive in a non-perturbative way equations for the effective permittivity 
tensor. The second approach is the result of an expansion in powers of puj Pi \pm (pi) is the impedance of the metallic 
(dielectric) component], combined with scaling ansatzes near pc. 

First of all, both approaches are consistent with each other, the only difference is that the scaling exponents 
predicted by the replica approach have their SEMA values. Second, the main result is the following: the scaling of 
the Faraday coefficient is the same for the weak field and the strong field regimes (as long as both uiej ^ ca/ and 
oj£iH ^ (Tj\/). In particular, we found that the Faraday angle (proportional to JF) is proportional to the magnetic 
field B even for strong field, as long as p < Pc and ^I'^^'i ^ 1- We can thus predict that it should be in principle 
possible to obtain large values of the rotation angle in such a system. For instance, for Ap of order 0.1 (which is 
realistic in experiments), £/ of order unity and (Tj\//w of order 100 for semiconductors, and with lj in the microwave 
region (the ratio (Tj\/ / ^ should not be too large since Eg is proportional to its inverse) , one obtains for the Faraday 
coefficient 

^ c lO-'H, (5.1) 

which can be made of order unity using currently available magnetic fields and high mobility doped semiconductors. 
One should recall that in homogeneous dielectrics, the Faraday coefficient is usually much less than 1: Typical values 
for a IT magnetic field, and for a wavelength in the visible spectrum (A ~ 0.6pm), are of order 10~^ for dielectrics 
like quartz, and of order lO^'^ for thin ferromagnetic metallic iron films. Measurements of the Faraday effect below 
Pc in a percolating metal-dielectric composite could provide an important test of the scaling predictions in both the 
strong field and the weak field regimes. Such experiments would have to involve either propagation or reflection of 
microwaves by a metal-dielectric composite with metallic inclusions that are smaller than the relevant skin depth. 
Both approaches also predict that the transverse diagonal elements of ie have an imaginary part that is proportional 
to Lue'j/aM, and are independent of H [up to terms of order 0{pm) in the high contrast expansion]. This means that 
there will be some dissipation. 
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APPENDIX A: 



We want here to solve the Maxwell equation satisfied by the electric field E (where — lu /c) 

[V X V X +e(r)fc^]E(r) = 0. (Al) 
The tensor e is a random variable equal to ii with probability p and to £2 with probability q = \ — p [ii and £2 



are tensors of the form given in (1.1)]. In order to obtain an integral equation for the electric field, we first write 
£(r) = £qI + (5£(r) (where £0 is an arbitrary constant which will disappear at the end of the calculation). One can 
then easily show that E(r) is also the solution of the following equation 

E(r) - Eo(r) + j dr'G(r - r')<5£(r')E(r'), (A2) 

where G is the dipolar tensor for the uniform medium of permittivity £0 and where the quantity Eq depends only on 
the boundary conditions and is assumed to be uniform. 
We will use the Fourier transform of Gapiv) given by 
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We work in the quasi-static limit which means that we take the hmit fco going to zero. In this case, the dipolar tensor 
is given by 

G.Mfc)=-^ for fc^O. (A4) 

£qK 

At fc = 0, the value of this tensor is Gq = —Sajj/ideo), where d is the space dimension, here equal to three. For 
values of fco that are too large, we cannot define an effective permittivity tensor and we have to introduce the notion 
of spatial dispersion (for a review, see e.g., Ref. |l6| and for a study using the replica method see Ref. [l^ ). 



Averaging equation (A2), after inverting it, and averaging it before inverting it leads to the following exact relation 

< Mj{k - 0) (1 - Gofe'e)-^, (A5) 

where is the inverse of the random operator 

M^0{r, r') = 6a,fi6{r - r') - (G(r - r')(5e(r'))a/3. (A6) 

Calculation of the effective permittivity tensor is thus reduced to a calculation of the average over the disorder of 
the inverse of a random operator, namely M . In order to do this, we will use the replica method which allows us to 
express the elements of (after using a Gaussian inversion formula) in terms of the functional integral 

M-^i(r,r') ^ / P(^,V.)C(r)^^(r')e/'^"-'^"-'^"^.»'^=^'-^^"^(^'^'^^^^^'\ (A7) 



where ip^ (with a = 1, n) (and its conjugate "0) are replicated Grassman fields satisfying the usual anticommutation 
relations 

{V(r), V(r')} = {V^(r), Vi(r')} = {^(r), ^^(r')} = 0. (A8) 



The limit n = is implicitely taken in (A7) and, as usual in the replica method, we first consider n as an integer 
and then take the limit n going to zero at the end of the calculation (without adressing the problem of analytic 
continuation) . 

It is now easy to average over the disorder, and we obtain 

{M-')Mr - r') = / V)C(r)V^?(r')e''% (A9) 
where the effective Hamiltonian is given by 

fie^ J drdr' ^ r^{r){6a^pS{r - r') - (G(r - r')<5£i)a;3}V^(r') 



a,/3,a 



drn In 



(AlO) 



The matrix A is equal to £2 — £1 and 77 = {1 — p)/p (note that rj is the inverse of A — p/q which appears in the 
main text). As usual in the replica method, the average over disorder introduces coupling between different replicas 
(if there is no coupling, then the |a¥|eraging is trivial) and in order to study this complicated effective Hamiltonian, 
we will use a variational principletSEJ. This principle consists of finding the best Gaussian approximation Tip to the 
effective Hamiltonian He- Denoting by K'^ the kernel of Ho (the variational approximation thus reads {M'^) ~ K), 
we have to minimize with respect to K the following variational free energy $ 

$(i^) =To + {He - no)o. (All) 

where Tq is the free energy associated with Ho and where (•)o denotes an average using Ho- We thus obtain the 
following equation 

{K-')o.p = 6c.p - {G{k)Sie)c.p, (A12) 
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with 



Se = ei+^{- 

m>l 

The tensor K can be inverted and we obtain 



K 



1 + ]S3 



k • q k (g) q 



D 



D 



(A13) 



(A14) 



where D — eofc^, q = <5e*k (^e* is the transpose of Sie) and where (g) denotes the usual dyadic product. We can then 
compute KG and we find that it is the dipolar tensor for the efTective medium 



«K/3 



It is then easy to integrate KG and we obtain 



dk 

(2^ 



{kG)^f-_ 



(A15) 



(A16) 



where /i is a diagonal matrix with diagonal elements fii = fJ.2 — l/se and /is = l/sze- The self-consistent equation 
( [A13| ) can thus be recast as the following matrix equation 



due 



{ee 



-uiil/3\ 



where the brackets still denote an average over the disorder. It is then easy to show that 

-«e2/3ee \ 





-I 


\ -He Se J 





due 



where £2 denotes the restriction of the tensor e to the {x, y) subspace 



£2 = 



e le 
-ie e 



(A17) 



(A18) 



(A19) 



The equation along the z— axis is decoupled from the preceeding one and is Eq. (2.1) of the main text. Using the 
following relation 



cosh ae —i sinh as 
i sinh ae cosh ae 



(A20) 



we obtain after simple manipulations equations (2.2) and ( |2.3| ) of the main text. It should be noted that all the 
calculation presented here can be used without any changes for either real or complex values of e and e. This justifies 
the analytical continuation of the formulas obtained in the framework of the Hall effect in order to describe the 
Faraday effect. 



APPENDIX B: 

The equations for the Green tensor &P^{r,r') can be solved in almost closed form if the system occupies all space, 
and if the resistivity tensor p is constant everywhere and its symmetric part is a scalar tensor, i.e., if 

p(r) • V = /5o • V = pqv + b X V (Bl) 

for any vector v. In that case, the Green tensor depends only on r — r', and we can define its Fourier transform by 



Gi'^\ci)^ J d(r-r')G[f;^(r-r')e-'>-C-'-'). 



(B2) 
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Using (3.3), it is easily found that this Fourier transform satisfies the following linear algebraic equations 

[poq^ - - pogo(q • + (b • q)(q X ^f)^ = 5^p, (B3) 

which can be solved to yield {£ai3-y is the basic antisymmetric tensor) 

n<^pa)f N _ (Pog^ - fc2)5a/3 - [poiPoq^ - k^) + (b • q)^] q^qfi/k'^ + (b • q)£a/3797 /T^.^ 
^"^^ - (po,2_fc2)2 + ^2(b.q)2 • (B4) 

Clearly, G^''^''(q) diverges in the limit fc ^ 0. However, if one calculates the Fourier transform of (V x G*'^°'')a, namely 

/ ^ Hp»)\ - (b ■ q)(g^ 5ap - q^qti) - [poq^ - k^)£ap^q^ , . 

^"^"^ (po<72-P)2+52(b.q)2 ' ^^^) 

then the limit fc ^ can be taken without any problems. That is why we had to include the term k"^ g''^^ in the 

equation for the Green tensor [see (3.3)], deferring the limit fc ^ until after the calculation of (V x G''^'')^. 

In order to investigate the symmetry properties of G*-''-' (r, r'), we use integration by parts or Green's theorem to 
get, for any vector fields A(r), B(r), and second rank tensor field /o(r), 

/ dr{A • V X [p • (V X B)] - B • V X [p • (V X A)]} 
Jv 

= - (f [{dS X A) • p • (V X B) - (dS X B) • /5 • (V X A)] 
Jav 

+ I dr [(V X A) • p • (V X B) - (V X B) • p • (V X A)] . (B6) 
Jv 

If p is symmetric, then the integrand in the last volume integral vanishes everywhere. Substituting 

^^(r)^G(f)(r,ri), S.(r) ^ Gg(r, r^) (B7) 



in this result, and assuming that p is a symmetric tensor, we thus get, using ( p.3[ ) and (3.4) 

Gi^j(ri,r2)=G(''j(r2,ri) if = p. (B8) 
Note that if p is non-symmetric, then in general G^''2(r, r') is not a symmetric kernel 
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